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Interaction of a domain wall with boundaries of a system is studied for a class of stochastic driven 
particle models. Reflection maps are introduced for the description of this process. We show that, 
generically, a domain wall reflects infinitely many times from the boundaries before a stationary 
state can be reached. This is in an evident contrast with one-species models where the stationary 
density is attained after just one reflection. 

I. INTRODUCTION 

There are many intrinsically nonequilibrium phenomena which can be observed already in simplest systems of 
driven diffusing particles, a recent review of which can be found e.g. in 0, 0. Phase transitions induced by 
spatial boundaries of a system is one of those phenomena which was studied in detail for models with one-species 
of particles [3- 0- 0] and for some multi-species models • The ability of a nonequlibrium system to "feel" the 
boundaries constitutes a key feature of driven systems: in fact, the boundaries dominate the bulk giving rise to the 
phase transitions. It is a flux which brings information from the boundaries to the bulk: in absence of a flux, the 
boundary conditions play only a marginal role as one knows from equilibrium statistical mechanics. 

A boundary problem in a one-dimensional driven diffusive system can be formulated as follows: the system is coupled 
at the ends to reservoirs of particles with fixed particle densities. A dynamics in the bulk and at the boundaries is 
defined via hopping rates, which are time- independent. In this case, after a certain transition period, the system 
will approach a stationary state, characteristics of which (the average flux, the density profile, the correlations) do 
not depend on time. It is of of interest how the system relax to the stationary state. Since a stationary state is 
independent on an initial state, one can choose an initial condition to be a domain wall interpolating between the 
boundaries. It has long been recognized that boundary-driven phase transitions are caused by the motion of a shock 
0, so that the above initial choice is a very natural one. A key issue is to understand how the domain wall interacts 
with the boundaries. For the reference model in the field, an Asymmetric Simple Exclusion Process (ASEP) with 
open boundaries, a stationary density is reached after just one reflection, or interaction with a boundary. 

We studied the reflection for the model with two conservation laws introduced in , and suprisingly observed 
infinitely many reflections occuring before a stationary state can be reached. More precisely, a domain wall after the 
first right reflection (i.e., the interaction with the right boundary) attains a different value of the density which is 
then changed after the first left reflection, etc.. This process continues iteratively. The bulk densities after multiple 
reflections asymptotically approach the stationary values. 

Description of iterative reflection maps constitutes the main subject of the present paper. 

Although the observed phenomenon of infinite reflections should be general, we use for the study a class of driven 
diffusive models having a simple stationary state on a ring for the sake of notational simplicity and to eliminate 
possible source of errors which may result from approximating a bulk flux, boundary conditions etc.. The findings 
are supported with Monte Carlo calculations and hydrodynamic limit analysis. 

The paper is organized as follows: In section[H]we define the model and its symmetries and describe the stationary 
state on a ring. Hydrodynamic limit equations arc given in section ITTT1 In section llVl reflection maps are introduced 
and infinite iterative sequences are analyzed. Details concerning boundary rates can be found in the Appendix. 

II. THE MODEL 

The model that we shall take as the example for our study is a particular case of a more general model introduced 
in 00, which can be viewed as a two- line generalization of ASEP. There are two parallel chains, chain A and chain B 
consisting of N sites each. Any site can be empty or be occupied by a particle. The particle at a site k can hop to its 
nearest right-hand site fc + 1 on the same chain provided it is vacant. The rate of hopping depends on the occupancy 
of the adjacent sites on the other chain, see Fig. 2] Hopping between the chains is forbidden. The hopping rates are 
fully symmetric with respect to an exchange between A and B. Correspondingly, there are 4 different rates a, (3, 7, e 
(see Fig. QJ. However if one chooses them to satisfy a + (3 = 2j = 2e, then all particle configurations occur with the 
same probability in a stationary state (see 9J). In this case, the rate j3 can be used to parametrize all the rates, 



a = l; T = e=(l + /3)/2. 



(1) 
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Stationary fluxes on a ring with the densities p A and p B of the particles on the A and B -chain respectively are given 
in the thermodynamic limit (see [9( ) by 

j A ( P A ,p B ) = p A (l-p A )(l + (0-l)p B ) 

j B (p A ,P B ) = P B (l-p B )(l + (l3-l)p A ) W 

For = 1 all the rates Q become equal, which corresponds to an absence of an interaction between the chains, and 
the system splits into asymmetric exclusion processes 0, fToL ITU . Another choice of the rates a — = 7 ^ e was 
considered in 

The rates are by the definition nonnegative, so is in the range < < 00. However it is sufficient to consider 
< < 1 because of a particle-hole symmetry. Indeed, an exchange of particles and holes plus the substitution 
— > 1//3 leaves the model invariant. 

In our context, we shall deal with a finite chains, where particles can be injected at the left end and be extracted at 
the right end of the system. We choose the boundary rates corresponding to effective stationary reservoirs of particles 
with fixed particle densities. Explicit expressions for the boundary rates are given, see the Appendix, in terms of the 
boundary densities. We denote the densities of ^4(i?)-particles in the left boundary reservoir as p A (p B ), and p^(pfj) 
for the right boundary reservoir. 



III. HYDRODYNAMIC LIMIT 



A naive continuum (Eulerian) limit of our stochastic dynamics on a lattice hk(t) —* p (x,t), rhk[t) — > p B (x,t) is a 
system of conservation laws 

dp Z ( X ,t) dj Z (p A ,p*) 9V y A U W 

0*" + to ^to^ Z = AB ' 



p z (0,t)=p z ; p z {N 1 t)=p z R 

where on the right-hand side of © a phenomenological vanishing viscosity term has been added. This simplest 
regularization term leads to the correct answer for the Riemann problem, as compared to the stochastic model, see 
but fails to describe a reflection from the boundaries (see Table 0. An adequate viscosity term is obtained by 
averaging of exact lattice continuity equations of the stochastic process 

J^fc - Sk-i - t (4) 

^-mk=3k-i-3k (5) 

for occupation number operators h]~ , rhk , and taking the continuous limit (hk) — > p A (x,t), {rhk+i) — * p B (x + a,t). For 
the case the flux operator reads (see 

Jk = "fc(l - n k+ i) (1 + 1 {rhk + ™>k+i) J , (6) 



and j B is obtained by an exchange n *-* rh in the above. We substitute JBJ into JU,©, average, factorize and Taylor 
expand the latter with respect to a lattice spacing a 1. In the first order of expansion, one obtains 

<l + V -l)f)°£) (7) 
+ (8) 

A numerical integration of Q , JSJl shows that the dynamics of the stochastic process is described adequately. At the 
same time, predictions of © strongly disagree with the stochastic dynamics (see e.g. Table[IJ|. Thus, in contradiction 
to the one-species case, where a choice of the vanishing viscosity is very robust, for a model with two and more species 
the specific choice of the viscosity becomes crucial: different choices give different answers. 
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IV. REFLECTION MAPS 



We explain our approach firstly for noninteracting chains (3 = 1. In this case th e dy namics of each chain reduces to 
the one of the ASEP with an injection rate pl and an extraction rate (1 — p R ) , see [l(J, El > corresponding to reservoirs 
of particles with densities pi, and pr on the left and on the right boundary, see the Appendix. To study interaction 
with the right boundary, choose an initial state to be a homogeneous state (fik) = Pbuik = Pl, for all fc, matching 
perfectly the left boundary. Consequently the particle distribution at the left boundary will not change in time, while 
at the right boundary there is a mismatch if p^ ^ pr, which has to be resolved. Monte-Carlo simulation shows that 
one of the following scenarios is realized, see Fig. [21 ( a ) a thin boundary layer develops interpolating between the 
bulk density pbuik = Pl and pr, and stays always attached to the right boundary, (b) a shock wave of the density pa 
develops and propagates to the left boundary, (c) a rarefaction wave with the density p = max(pR, 1/2) forms and 
spreads towards the left boundary. 

Analogously, we study the left reflection, choosing an initial condition p^ ^ pbuik — PR- Then, the mismatch will 
be only at the left boundary. We can predict interaction results by exploiting a particle-hole symmetry of the model: 
an exchange of particles and holes and the left with the right boundary leaves the system invariant. Hence a dynamics 
of the model with p^ = Pb u ik PR is equivalent to a dynamics of the dual model 

Pl + Pbuik = Pr, wh ere p' L = 1 - Pr, p' R = 1- p L : (10) 

The corresponding local densities of the dual model satisfy 

(n k ) = {n'N-k+i); k = l,2,...JV. (11) 

As an example, Fig. |2i shows the Monte Carlo evolution corresponding to shock propagation, dual to the one on 
Fig. 03, through the transformations (|10|) . . The other scenarios for the left reflection can be deduced from 
Fig.|2i,c, us i n g (dnj , ((TU) . 

We shall classify an outcome of a reflection from a boundary by an average particle density p re fi of the reflected 
wave measured at some point k* not too close to the boundary to avoid boundary effects. The density p re fi is 
measured after the reflection wave has passed the point k* but before it has reached the other boundary, to avoid 
possible interference. For instance, reflected wave has the density p re fi = PL for the scenario (a), prefi = PR for the 
scenario (b) and p re fi = 1/2 for the scenario (c) on Fig. El 

We can summarize the results of all possible reflections from the right boundary pr by plotting the density of the 
reflected wave versus the initial bulk density, see Fig. [3] We shall call this type of graph a reflection map. Reflection 
map for the left boundary is obtained using i|l(J|) . l|ll[l . Comparing the reflection maps with the stationary phase 

diagram of ASEP dsim, 

one can make an important observation. Namely, the stationary state density is achieved 
after one reflection from any of the boundaries. We have checked this statement to be true also for a two-parameter 
model with a next nearest neighbour interaction (KLS model), the stationary phase diagram of which is obtained in 
1 

In the following we show that for a two-component system it takes infinite number of reflections until the stationary 
density is reached. 



A. Two interacting chains 

Consider the two-chain model JQ) coupled to boundary reservoirs with the densities of A- and B- particles 
and p R , pjj at the left and at the right end respectively. Consider a reflection problem by choosing an initially 
homogeneous particle distribution with the densities 

Pbuik = Ph Pbuik = Pl for the right reflection (12) 



Pbuik = Pr\ Pbuik = Pr for the left reflection (13) 

We shall classify results of a reflection as it was done above for ASEP. Denote the densities of the particles in the 
reflected wave in A— and _B-chain as r A ,r B . Suppose that after interaction with a boundary an interface appears 
between Pbuik' Pbuik ano - rA > rB moving towards the other boundary. Due to particle number conservation in the bulk, 
a Z-componcnt of the interface moves with the velocity |13| 

AZ _ aZ 

yZ = Jbulk Jr ^ , 14 j 
Pbuik " r 
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where we used the short notations j z nlk — j Z (p A u \ k , Pbuik)> 3x ~ J Z { rA i rB )- Since there is an interaction between 
the chains, the velocities must coincide in both chains because a perturbation in one chain causes the response in the 
other and vice versa. This gives a restriction V A = V B , or 

aA _ aA -b _ -b 

Jbulk Jr _ Jbulk -h (15) 
Pbulk ' Pbulk ' 

defining implicitly the allowed location of the points r A ,r B . 

For demonstration purposes, we choose j3 = in Q which corresponds to the maximal interchain interaction. This 
choice simplifies drastically analytic expressions and at the same time preserves qualitative features of the general 
case. For (3 = 0, Ea. l|15() has two families of solutions, given by 

r A /r B = P U/pU (16) 
(1 - r A )(l r B ) = (1 - ^ ulk )(l - p B ulk ) (17) 

It turns out that the densities of waves resulting from the interaction with the right boundary satisfy {T5|| while those 
resulting from the left reflection satisfy (|T7jl . 



1. Right reflection 

Consider the right reflection l|12|l first. Monte Carlo simulations, independently confirmed by a continuous model 
integration, lead to the following conclusion: 

the reflection densities r A , r B either coincide with the p A uik , Ptmik ( m this case boundary layers analogous to Fig.Et 

develop), or depend only on the ratio 7 = b w ■ In the latter case, r B ("f)/r A (j) = 7, sec (1161) . On Fig. 01 a typical 

Pbulk 

dependence of r A on the initial bulk density along the line = 7 is shown. Analogously to the left graph on Fig. 

fbulk 

a discontinues change of r A occurs, defining a point of a first order boundary driven phase transition |l4l |. 

Collection of points r A (7) , t b (7) for all possible initial bulk densities < 7 = < 00, constitutes a continuous 

^bulk 

curve TZ shown on Fig. |SJ The shape of TZ depends only on the right boundary densities p^, pj| . In the following we 
derive some analytical properties of TZ, namely the coordinates of three points on TZ and the derivative in one of them. 

Trivially TZ contains the point Pr, Pr corresponding to the perfect match with the right boundary. Then, the 
location of end points of TZ can also be derived. Indeed, consider the system with the empty i?-chain: p B uik = 0. 
Then, the model becomes effectively an ASEP, with the rate of extraction of the particles i° xt given by (A. 9), which 
corresponds to the effective right boundary density in ASEP pn = 1 — t® xt . Consulting the reflection graph for ASEP 
Fig. |3 we find that corresponding reflected wave in the A-channel has the density max(l — i° xt , 0.5). Thus, the curve 
TZ has a point with the coordinates (max(l — (1 — /3^)(1 — ^p^), 0.5), 0). Repeating the arguments for the empty 
A-chain p A ulk = 0, we find the other end point on TZ, (0, max(l — (1 — p B )(1 — \pfy, 0.5)). 

Now, the exact derivative of TZ at the point p^, p 3 , can also be derived. Indeed, it is shown in the section ITTT1 that 
our system is described in the continuous limit by the system of equations J7JSJ). After the initial period of strong 
interaction with the right boundary, the distribution of particles near the boundary does not change anymore (see 
e.g. Fig|3b), and can be described by a time-independent solution of I|7I8|I . The current in the region will then be 
equal to that of the reflected wave. Integrating the time-independent part of Eq.{7|) for (3 = 0, one obtains: 

j A (p A ( X ),p B ( X )) ^l t (r A ,r B ) (l-p B (x)) (18) 

where the current of the reflected wave i^t a t(^ A , ^ B ) is the constant of integration. 

Suppose the result of the reflection (r A ,r B ) differs only infinitesimally from the right boundary values, r z w p^, 
Z = A, B. The Taylor expansion of Ijl8|l around the r A , r B gives in the first order approximation 

^)/ / + (^)/' B = A(i/<w ° )) (19) 

where the subscript R denotes evaluation at the point p^,p B , and we suppose the profile p z {x) to approach expo- 
nentially the values r A ,r B from the right. An equation for the other flux component j B , obtained analogously, can 
be combined together with l|19|) into a modified eigen-value equation 

8p B ) a \5p B )' * { l-p£ 
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(A B\ 

(Pr Pr) 


Theor 


Numerical 


"Naive" viscosity choice 


(0.2 0.8) 


-0.151388 


-0.15(5) 


-0.25 


(0.3 0.8) 


-0.190476 


-0.19(1) 


-0.285714 


(0.5 0.5) 


-1 


-1 


-1 


(0.4 0.6) 


-0.548584 


-0.55(3) 


-0.666667 


(0. 0.6) 


-0.2 


-0.2 


-0.4 


(0.6 0.) 


-5 


-5 


-2.5 


(p A , 1) 












Table I: Comparison of the analytical formula 121H for derivative of 1Z at the point (p'r, Pr) with the numerical results. 



where Vj is the Jacobian of the stationary flux © with the elements T>uy = dj u /dp v . For any given value of the 
boundary densities, A is determined self-consistently, and the ratio 8p B /8p A is the derivative of the curve 1Z at the 
point p^, p 3 ,. The system (I20II has two solutions, one is relevant for the right boundary reflection and one for the left 
reflection. The solution, relevant for the right reflection, reads for f3 = 0: 

Sp B (i - p b b)(p b r p a r - y^) 2 + K) 2 ~ p a rP b r) ^ 
Sp A 

We observe an excellent agreement between the theoretical and numerically evaluated values of j^-k, see Table [|] The 
last column of the Table gives the value of jp$- which would follow from the "naive" viscosity term as written in (J2J), 
which leads correspondingly to B = I in 11201 . and to the expression (t^x) = — \ p " Hoi . One sees from the 

v op 'naive ^~Pr 

Table that the latter choice gives a wrong prediction. 



P&0--P&) 



2. Left reflection 

Consider the reflection at the left boundary (|13|) . We observe that the densities of the reflected wave r A , r B either 
coincide with the left boundary densities p A ,p B or depend only on the value T characterising the respective solution 

(i-Pbuik)(i-Pbuik) = r (22) 

More precisely, the curves < T < 1 , see Fig. span the whole parameter space < Pb u ik! Ptmik — ^ 

into three regions & S oiid,$ dotted (containing solid and dotted part of curve C respectively), and & empty, containing 
the resting bottom left part of Fig. |5| In the region <E> empty, the wave matching perfectly left boundary reservoir, 
(r A , r B ) = (p A , p B ), is produced. In & S oHd, reflection wave densities lie on C and for any T are given by an intersection 
of C with the hyperbola (|22|l . Note that this intersection point is unique for any T. Finally, in Q dotted, the reflection 
wave consists of two consecutive plateaus. First plateau from the right has the density X A 7 X B (point of intersection 
between C with JUJ), and the second matches the left boundary. Since the velocity Vlx of the interface between the 
two is positive for all X A 1 X B belonging to the dotted part of reflection curve £, 

3 A (X A ,X B )-j A (p A ,p B ) j B (X A ,X B )-j B (p A ,pE) ^ n 
VLX = X^-pJ = XB-^p* > °' (23) 

the effective reflection result is (r A ,r B ) = {p A ,p B ) like in the region & empty At the point K where the dotted 
and solid part of the curve C meet, the velocity (|23|) vanishes, Vlx(K) = 0. Inside the region & S oUd the velocity 
Vlx < meaning that the intermediate shock interface stays "glued" to the left boundary and only one reflection 
wave survives. Note also that due to (J2HJ), the points X A , X B from the dotted part of C satisfy X A / X B — p A /p B as 
follows from fity . fity . 

Generically, the curve £, corresponding to arbitrary boundary densities p A ,p B , contains this point itself, corre- 
sponding to a perfect match with the left boundary. Secondly, it ends at the point (1, 1) for the case /? = which we 
consider. Indeed, if p^uik ~ 1> ^ ne cham B is almost completely full, so that no flux can flow at the chain B because 
of exclusion rule. But, a flux is very small also in the chain A, because the dominating configurations, see Fig^ 
have zero bulk hopping rate (3 = 0. On the other hand, the particles can enter with a finite rate (|A.3p . In the limit 
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lk — > 1 this yields the reflected densities (r A ,r B ) - (1, 1). Or, in terms of the reflection map, (Pb ulk , 1) — > (1, 1) 
Exchanging the role of the A- and S-chain leads to the same result, (l,Pbuik) ~~ * (•"■> 1)> due t° complete symmetricity 
of the rates. 

The reflection maps Figs l5l6l can be interpreted also in the following sense: the right reflection map Fig. (SJ shows 
the location of all possible stationary densities of a system with the given right boundary densities Pr,Pr- Indeed, 
after a right reflection the density of particles in the system is located either on 1Z or in the bottom left (low density) 
region. Correspondingly, a stationary state (whatever are the conditions on the left boundary) must belong to the 
same domain. Analogously, the left reflection map Fig. El considered alone, shows that the stationary densities of the 
system with the given left boundary densities either lie on the curve £ or match the left boundary p^jPl- 

We shall not explore further the analytical properties of the curves C and 7Z, but investigate what happens if we 
fix the right and the left boundary corresponding to the reflection maps Fig. El and Fig. El Initially empty system 
fills with the particles from the left end, with the densities, fitting the density of the left reservoir, Fig. EI Thus the 
wave Pb u ik i Pbuik = PltPl starts to propagate in the system towards the right boundary. Reaching the boundary, a 
reflected wave forms, with the densities r A given by the intersection of the line y/x = p£ /p^ with the curve 
1Z of Fig. El This reflected wave hits the left boundary then, resulting in the next reflected wave with the densities 
rf ,rf . The left reflection is controlled by the curve L at Fig. El therefore the densities of the reflected wave rf 
,rf are given by the intersection of hyperbola (1 — y)(l — x) = (1 — rjf)(l — rjf ) with the curve C. Since the curves 
1Z,C do not coincide with any of the curves defined by l|17llt)l) in our example, this process continues forever, though 
converging to stationary state S, the point of the intersection of C with 1Z, see Fig[7| After 2/c — 1 reflections, the 
wave with the densities i r 2k-i wm n ^ the right boundary, producing the reflected wave of the densities 



2k — 1 ' 2k— 1 VV111 UXXC UUUllUCllJ, ^lUUU^lllg ICllC^tCU VVCLVC Ul UllO u^noiuico '2k 

> r 2k> corresponding to intersection of the line y/x — i"2k-\l r 2k-\ with the curve 1Z. The latter wave, hitting the left 
boundary, produces the next reflected wave with the densities r^ fc+1 , the point of the intersection of hyperbola 

(1 — y)(l — x) = (1 — 7'^ £ )(1 — r^,) with C. This process is depicted schematically on FigJH] Several remarks are in 
order. 

• The densities {r^ ,rjf })£L of the reflected waves constitute converging sequence. Convergence is exponential, 
8u n +2k = e~ Kk 5u n , for n — > oo where we denote by 8u n the deviation from the stationary density at the n-th step. 
e~ K < 1 is a constant depending on tangential derivatives a, ft, <x\,j3i of the curves 1Z,£ and the characteristic 
curves, respectively, at the stationary point S: e"« = (1 - ^)(1 - ^)/((l + W^ 1 " N ° te 
that if at least one of the characteristic derivatives happens to coincide with a%, /3\ : a — a\ or /3 = f}\ then the 
sequence converges in one step. 

• The velocities of the reflected wave for big k converge to the finite characteristic velocities computed at the 
stationary point 5, as eigenvalues of the flux Jacobian, see 0. 

• The stationary densities are not reached at any finite time [l6|. 

We would like to stress that convergence to a stationary state through infinite number of reflections is a generic 
feature of a two-component model. However, in specific cases the convergence to stationary state can be achieved 
after finite number of reflections. E.g., if the left boundary densities are equal, the curve C will be a straight line y = x 
coinciding with the one of the solutions l|16(l , and consequently the corresponding stationary state will be reached in 
one step, like in the one-species case. 



V. CONCLUSION 



To conclude, we have studied two lane particle exclusion processes on the microscopic level, focusing on interaction 
of domain walls with the boundaries of the system. It was shown that a shock front interacts infinite number of 
times with the boundaries of the system in order to establish a stationary state density. It was also shown that 
hydrodynamic limit with a usually taken approximation of diagonal vanishing viscosity ( see, e.g., [TJ) fails to 
describe the microscopic dynamics. An adequate hydrodynamic limit is constructed. 

For the description of the reflection of domain walls from the boundaries, we introduced reflection maps, which 
show the locus of the allowed densities of the reflected waves. Careful analysis of such maps is necessary in order to 
solve a boundary problem for models with more than one species of particles. 
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Appendix: BOUNDARY RATES 

Here we consider our model on two parallel chains of length N where constant densities of particles are kept on 
the left and on the right boundary. This can be achieved by choosing appropriately boundary rates for injection and 
extraction of the particles. 

The average current through the left boundary of chain A is by definition 



where r. 1 ^ (A) is the probability to inject a particle on the first site k = 1 of chain A provided adjacent site is occupied 

(empty). On the other hand, the average current through the link between the sites k and k + 1 of chain A is given 
by 

j A = l{n k (l - m k )(l - % + i)(l - m k+1 )) + /3(n k m k (l - n k+1 )rh k+1 ) 

+ — ^ 1 (h k m k (l - %+i)(l - fhk+i)) + 1" - m k )(l - h k+1 )fh k+1 ) (A. 2) 

Let the left boundary density to be p A (/?£) for chains A (B) respectively and assume the same stationary density in 
the bulk. Then the bulk current and the boundary current must be equal j A = j A . Comparing IjA.ljl and (| A. 2|> . and 
accounting for the absence of correlations in the stationary state, see (9J, we obtain the injection rates 



r^(A) = <n)(l - m) + ^yl<n)(m> = p A (l + L-lrf )\ 



(A.4) 



Now, due to the symmetricity of the model, the injection rates r^j (B) ,rP n j (B) for the chain B are obtained by an 
exchange of p A <-> p£ in the above expressions: 

4j( iJ )=P£( 1 + ^ + (A-5) 
rL( B )=P B L (l + ^ 1 pi)) (A.6) 



1Iy v , ri, ^- 2 

Analogously, the current through the right boundary is due to the extraction of the particles from the right end k = N 
of the chains A and B, 

3R = ^rtC^XMl - m N )) + tl xt (A)(h N rh N }, (A.7) 

i| = t° xt (S)(m JV (l - njv)} + tg X i(B)(n N rh N ). (A.8) 
Assuming the right boundary densities (and the bulk densities) to be p^, p^ and equating l|A.7|) with 1A.2|) . we obtain 

i° xt (A) = /3(1 - - m) + ^±1(1 - n)(m) = (1 - p A ) (l + ^pf) (A.9) 
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^ „B 



4 xt (A) = (1 - pi) ( 1 + ^(1 + pg) ) . (A.10) 



The expressions for the extraction rates from the chain B are obtained by exchanging <-> p^ in the above, 

i° xt (B) = (l-pf)fl + ^-lp^ (A.ll) 



I) (1 + ^4 

tl xt (B) = (1 - p*) (l + ^1(1 + pi) ) . (A.12) 



For any fixed choice of boundary densities the model will relax to unique stationary state, independent on the 
initial conditions. For equal left and right boundary densities, pi — p A — p A , p 3 , = pf = p B , our definition of the 
boundary rates guarantees that the stationary state will be the uncorrelated homogenuous ( product) state (h^) = p A , 
(rhk) = p B for any number of sites N > 1. 
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Here we assume that the scenario of the type depicted on Fig. [5] 6, d takes place. 

With a different choice of p^,p|j, a reflection graph analogous to the right graph on Fig. EJcan be obtained. 

[15] If B = /, 120H becomes an eigenvalue equation for the Jacobian (Dj) Corresponding value of is then the derivative 

along the relevant characteristic at the point (pr,Pr). 
[16] Strictly speaking, the stationary density is reached exponentially in time with characteristic time proportional to the length 
of the system N. On the contrary, in one-species models like ASEP, the corresponding characteristic time is of order of 
l/(rate of hopping). 



Figure 1: The four elementary hopping processes, shown here for one of the chains, and their rates. In the study, the rates 
are chosen to satisfy 
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Figure 2: Scenarios of a domain wall reflection from a boundary, as given by Monte Carlo simulations. In all the cases, an 
initial distribution was constant. The average profiles after 150 and 300 Monte Carlo steps are shown. Averaging over 5 * 10 
histories is made. The graphs (a)-(c) show the reflection from the right boundary. The graph (d) demonstrates the reflection 
from the left boundary dual in the sense of llOH . lllll to the graph (b) 




Figure 3: The right boundary reflection map for ASEP 
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Figure 4: Right boundary reflection map for the two-chain model along the line Pbuik/Pbuik = 7- The boundary densities: 
pu = 0.2, pH = 0.8. Point 0(7) denotes the point of phase transition, where the bulk and the reflected fluxes become equal 
j z (a, r ya) = j z (rf ,77V). The locus of all points (r^,7r» ) and all points 0(7), 70(7) constitutes the upper and the bottom 
curve at Fig. |5J respectively. 



Figure 5: Complete right boundary reflection map for the two-chain model, showing the locus of the allowed densities of the 
waves reflected from the right boundary. The right boundary densities are: p^ = 0.2, p^ — 0.8. Squares show the Monte-Carlo 
calculation points, the upper continuous curve 1Z is obtained by integrating the hydrodynamic equation 17181 . The crossmark 
shows the right boundary densities. The bottom left (LD) region marks the low density phase and the bottom curve marks the 
line of the first order phase transition, see also Fig. 2] 
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Figure 6: Complete left boundary reflection map for the two-chain model. The left boundary densities are: p^J = 0.2, p£ = 0.4. 
The bottom left (LD) region marks the locus of initial densities resulting in the reflection wave which matches the left boundary 
(point L ). Squares show the Monte-Carlo calculation points. Dotted curves show hyperbolas 11221 for different V. 




Figure 7: Curves 1Z and £ from the Figs. I5lb1 combined together. Their crossing, S, is the stationary state density achieved 
through the infinite number of reflections. The crossmarks near the symbols 1Z and £ indicate the left and right boundary 
densities, and the star at the crossing — the stationary density obtained by Monte Carlo simulation of the system of 300 sites. 
The systems was equilibrated for = 4 * 10 5 Monte Carlo Steps (MCS), after which the averaging over MCS = 4 * 10 5 and 10 
different histories was done. 



12 




Figure 8: Approach to stationary point 5 through an infinite sequence of reflections. The numerated filled circles show the 
location of subsequent domain wall densities r A ,r B . The circles k, k + 2, k + 4 lie on the curve £ and correspond to the result 
of the left reflection, while those lying on the upper curve 1Z, correspond to the right reflection 



